ENERGY-EFFICIENT JOINT ESTIMATION IN SENSOR NETWORKS: 

ANALOG VS. DIGITAL 
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ABSTRACT 

Sensor networks in which energy is a limited resource so that 
energy consumption must be minimized for the intended applica- 
tion are considered. In this context, an energy-efficient method for 
the joint estimation of an unknown analog source under a given 
distortion constraint is proposed. The approach is purely analog, 
in which each sensor simply amplifies and forwards the noise- 
corrupted analog observation to the fusion center for joint esti- 
mation. The total transmission power across all the sensor nodes 
is minimized while satisfying a distortion requirement on the joint 
estimate. The energy efficiency of this analog approach is com- 
pared with previously proposed digital approaches with and with- 
out coding. It is shown in our simulation that the analog approach 
is more energy-efficient than the digital system without coding, 
and in some cases outperforms the digital system with optimal 
coding. 

1. INTRODUCTION 

A typical Wireless Sensor Network (WSN), as shown in Fig. Q 
consists of a fusion center and a large number of geographically 
distributed sensors. The sensors typically have limited energy re- 
sources and communication capability. Each sensor in the network 
makes an observation of the target of interest, generates a local sig- 
nal (either analog or digital), and then sends it to the fusion center 
where the received sensor signals are combined to produce a final 
estimate of the observed signal. Sensor networks of this type are 
well-suited for situation awareness applications such as environ- 
mental monitoring and smart factory instrumentation. 




Fig. 1. Sensor network with a fusion center. 

Decentralized estimation has been studied first in the context 
of distributed control Q, distributed tracking 0, and most re- 
cently in wireless sensor networks O. Among these studies, it is 
usually assumed that the joint distribution of the sensor observa- 
tions is known. In practical systems, the probability density func- 
tion (pdf) of the observation noise is hard to characterize, espe- 
cially for a large scale sensor network. This motivates us to devise 
signal processing algorithms that do not require the knowledge of 
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the sensor noise pdf. Recently, universal decentralized estimation 
schemes (DES) without the knowledge of noise distribution have 
been proposed in jvl and |8l. In (tJ, the author considered the uni- 
versal DES in a homogeneous sensor network where sensors have 
observations of the same quality, while in (8|, the universal DES 
in an inhomogeneous sensing environment was considered. These 
proposed DESs require each sensor to send to the fusion center 
a short discrete message with length decided by the local Signal 
to Noise Ratio (SNR), while the performance is guaranteed to be 
within a constant factor of that achieved by the Best Linear Unbi- 
ased Estimator (BLUE). An assumption in these proposed schemes 
is that the channels between sensors and the fusion center are per- 
fect, and all messages are received by the fusion center without any 
distortion. However, due to power limitations and channel noise, 
the signal sent by each individual sensor to the fusion center will 
be corrupted. Therefore, the goal of the transmission system de- 
sign for the joint estimation problem is to minimize the effect of 
the channel corruption while consuming the minimum amount of 
power at each node. 

Historically, if the sensor observation is in analog form, we 
have two main options to transmit the observation from the sen- 
sors to the fusion center: analog or digital communication. For the 
analog approach, we keep the observation signal analog and fur- 
ther use analog modulation schemes to transmit the signal, which 
is also called an amplify-and-forward approach. In the digital ap- 
proach, we digitize the observation into bits, possibly apply chan- 
nel coding, then use digital modulation schemes to transmit the 
data. It is well known (IH, O) that for a single Gaussian source 
with an AWGN channel, the amplify-and-forward approach is op- 
timal. However, for an arbitrary source with multiple observations 
and multiple transmission channels, it is unclear which approach 
will lead to a smaller power consumption while meeting the dis- 
tortion requirements at the fusion center. 

Another severe challenge facing sensor networks is the hard 
energy constraint. Since each sensor is equipped with only a small- 
size battery, for which replacement is very expensive if not impos- 
sible, the power consumption must be minimized to increase the 
network lifetime (6). Therefore, energy efficiency can be used as a 
performance criterion to evaluate different sensor network designs 
under the same distortion requirement. In (9|, a digital system is 
proposed to minimize the total power consumption for the joint es- 
timation problem, where each sensor quantizes the analog observa- 
tion into digital bits and transmits the bits to the fusion center with 
uncoded MQAM. The number of quantization bits for each sensor 
is optimized with the target of minimizing the total transmission 
power across all the sensor nodes to achieve a given distortion. In 
this paper, we propose an analog counterpart for the same system 
and compare the energy efficiency between the analog approach 
and the digital one. As in 0, we assume that the observed signal 



is analog and bounded, the fusion center deploys the best unbiased 
linear estimator, the observation noise is uncorrelated across dif- 
ferent sensors, and only the variance of the observation noise is 
known. 

Our paper is organized as follows. Section II discusses the 
problem formulation. Section III compares the energy efficiency 
between the analog approach and the digital approach via some 
numerical examples. Section IV summaries our conclusions. 

2. MINIMUM POWER ANALOG INFORMATION 
COLLECTION 

We assume that there are K sensors and the observation Xhit) at 
sensor k is represented as a random signal 6{t) con'upted with the 
observation noise nk{t): Xk(t) = 6(t)-\-nk{t). Each sensor trans- 
mits the signal Xk{t) to the fusion center where 9{t) is estimated 
from the Xfc(t)'s, k = 1, - ■ ■ ,K. We further assume that nk{t) is 
of unknown statistics and the amplitude of x{t) is bounded within 
[—W, W], which is defined by the sensing range of each sensor. 
For simplicity we assume W = 1, but our analysis can be easily 
extended to any values. We also assume that the network is syn- 
chronized, which may be enabled by utilizing beacon signals in a 
separate control channel. 

We assume that Xk{t) is also band-limited and the informa- 
tion is contained within the frequency range [—B/2,B/2]. We 
consider an analog Single Side-Band (SSB) system with a coher- 
ent receiver II II . The transmitted signal is given by 

yt{t) — 2\/q cos {LJct)x{t) + 2\/Qsin {LJct)x{t), (1) 



for which the average transmission power is 

P = 4aP^ < 4aW^ 



(2) 



where ia is the transmitter power gain and Px is the peak power 
ofx{t). 

The received signal at the fusion center is given by 

yr{t) = COS {ujct)x{t) + 2^ygsin {ujct)x{t) + nc{t) 

(3) 

where g is the channel power gain, x{t) is the Hilbert transform of 
x{t), and nc{t) is the channel AWGN. Hence, at the output of the 
coherent detector the signal is 

y{t) = ^^x{t) + -ni{t) cos{Tr—t) + -nf (t) sm(TT — t), 

(4) 

where ni{t) + jn^{t) is the complex envelope of ndt). After 
passing y{t) through a low-pass filter and sampling the baseband 
signal at a sampling rate B, we can obtain an equivalent discrete- 
time system. Since we have K such sensors, the overall system 
is shown in Fig.|2| The K transmitters share the channel via Fre- 
quency Division Multiple Access (FDMA), which has the same 
spectral efficiency as the Time Division Multiple Access (TDMA) 
that is used in (9|. 
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The received signal vector at each time instance is given by 

y = he + V, (5) 
where . , + 

y = [yi,j/2,--- ,t/x]' , 

V = [yaigmi + rici, • • • , y/OLKgRUk + Uck]^ , 

and f means transpose. 

According to IIOI . the best linear unbiased estimator (BLUE) 
for 9 is given by 

,tT3-ll,l-lutT3-l 
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where the noise variance matrix R is a diagonal matrix with Rkk = 
a\a.kgk + ^fe with cr^ the variance of the sensor observation noise 
nk{t), k = 1, • ■ ■ , K. The channel noise variance ^1 is defined 
by the noise power spectral density and the bandwidth B. 
The mean squared error of this estimator is given as IIOI 
Var[6] = [h'l'R^^h]"^ 
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(7) 



According to Eq. J2j, the transmit power for node k is bounded 
by 4W^ak- Therefore, the minimum power analog information 
collection problem can be cast as 

K 
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where Dq is the distortion target. However, this problem is not 
convex over the 's. 
Let us define 

1 
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Then the above optimization problem is equivalent to 
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where we see that the variable Ok can be completely replaced by a 
function of Vk- Therefore, the problem can be transformed into a 
problem with variables {ri, r2, . . . , ta-} shown as follows: 
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(8) 



Fig. 2. Amplify and Forward 



which is convex over ru - The upper limit on in the second con- 
straint is due to the fact that = 1 / ( a? + ) and > 0. 
Now we solve Eq. Js}. Its Lagrangian G is given as 
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(9) 



for < Tk < —J- , Vfc, which leads to the following Karush-Kuhn- 
Tucker (KKT) conditions 1121 : 



9k (1 



K 
k = l 



Ao = 0, 



Vfc 



Do 



= 



for < Tfc < -^,\/k. Without loss of generality, we rank the 
channel quality such that 



< ^ < . . . < ^ where we call the 

ffl — 92 — — SK 

quantity fj- the channel SNR, and we define 



^A(M) 
B{M) ' 



for 1< M < JC, 



(10) 



where A{M) = ^ -^-^ and B{M) 



M 
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Let us find Kx such that /(ifi) < 1 and /(isTi + 1) > 1. 
Using the same techniques as in (9), we can show that this K\ is 
unique unless f{M) < 1 for all 1 < M < K, in which case we 



take Ki 
and 



K. Then the KKT conditions give Ao 



WA{Ki) 



opt 




Vfc 



(11) 



where (x)^ equals when x < 0, and otherwise is equal to x. 
Hence, by definition, we have 



opt 



opt 



5^ 1 - ^^r* 



,Ki, (12) 



and Qfc = otherwise, where rjo = A{Ki) / B{Ki). 

Therefore, the optimal power allocation strategy is divided 
into two steps. In the first step, a threshold for k is obtained accord- 
ing to Eq. <10t . For channels with SNR worse than this threshold, 
the corresponding sensor is shut off and no power is wasted. For 
the remaining active sensors, power should be assigned accord- 
ing to Eq. <12t . From Eq. <12> we see that when the channel is 



fairly good, i.e.. 



tVo ^ li we have a?^ oc 



which 



means the optimal solution is inversely proportional to the square 
root of the channel SNR. When ^|t'7o is close to one, the opti- 
mal solution may no longer have such properties. For all channel 



conditions, the power is scaled by the factor -ij, which means that 

more power is used to transmit the signals from the sensors with 
better observation quality. 

We now solve the optimization problem for some specific ex- 
amples. We assume that the channel power gain ~ where 

dfc is the transmission distance from sensor k to the fusion center 
and Go ~ —30 dB is the gain at d = 1 m. As in (9), we generate 
al uniformly within the range [0.01, 0.08]. We take B = 10 KHz 
and = —90 dBm, fc = 1, • • • ,K. For an example with 100 sen- 
sors, Fig.|3|(a) shows the relative power savings compared with the 
uniform transmission strategy where all the sensors use the same 
transmission power to achieve the given distortion target. The rel- 



\/Var(d) 
E(d) 



ative power savings is plotted as a function of 7? = 
the distance deviation normalized by the mean distance. For each 
value of R, we average the relative power savings over 100 random 
runs where in each run the dk's are randomly generated according 
to the given R. As expected, a larger variation of distance and cor- 
responding channel quality leads to a higher power savings when 
this variation is exploited with optimal power allocation. 

In Fig.|3|(b), the number of active sensors over R is shown for 
an example with 10 sensors. We see that when the transmission 
distances for different sensors span a wide range of values (i.e., R 
is large), more sensors can be shut off to save energy, since the 
remaining sensors have very good channels. 




(a) Normalized distance deviation (R) 



(b) Normalized distance deviation (R) 



Fig. 3. (a) Power savings of optimal power allocation vs. uniform 
power; (b) Number of active sensors versus distance deviation 

In our model we minimize the power sum P^, i.e., the L^- 
norm of the transmission power vector P = (Pi, P2, . . . , Pk)- 
If the channel gain and the variance of the observation noise for 
each sensor are ergodically time- varying on a block-by-block ba- 
sis, minimizing the L^-norm of P in each time block minimizes 
Pk} with E{} the expectation operation. In other words, 
it maximizes the average node lifetime, which is defined as 
'k'^k Ei'p^.} ^^^^ Eo the battery energy available to each sensor 
(we assume that Eq is the same for all the sensors). This can be 
proved by the fact that ^ Y.k "Bptl - e{^y: ^Pk} ' 
when the channel is static and the variance of the observation noise 
is time-invariant, minimizing the i^-norm may lead some individ- 
ual sensors to consume too much power and die out quickly. In 
this case minimizing the L°°-norm, i.e., minimizing the maximum 
of the individual power values, is the most fair for all sensors, but 
the total power consumption can be high. As in (9), we can make 
a compromise to minimize the L^-norm of P. In this way, we 
can penalize the large terms in the power vector while still keep- 
ing the total power consumption reasonably low. For the L^-norm 
minimization, the problem formulation becomes 



min 



s. t. 




which we can solve using interior point methods 1121 . In the next 
section, we compare the power efficiency of the analog and the 
digital approaches previously discussed, where we minimize the 
L^-norm of P. We use the L^-norm since it simplifies the com- 
parison, but the comparison can be made for any power norms. 

3. ANALOG VS. DIGITAL 

In order to transmit the observed analog signal with frequency 
range [—B/2, B/2] to the fusion center, each sensor in the dig- 
ital system proposed in (9) first samples the signal at a sampling 
rate B, then quantizes each sample into bk bits, and finally uses 
uncoded MQAM to transmit the bk bits with a symbol rate B 
and constellation size M = 2'''= . Therefore, the total transmis- 
sion bandwidth in the passband is approximately equal to KB for 
the digital system where TDMA is used for the multiple access. 
For the SSB scheme used in our analog system, each sensor only 
occupies B/2 in the passband such that the total bandwidth re- 
quirement is KB/2 when FDMA is used for the multiple access. 
Therefore, under the assumption of orthogonal channel usage, the 
analog system can support double the number of sensors that the 
digital system supports in the same amount of bandwidth. This is 
mainly caused by the fact that the digital approach proposed in (9) 
forces the transmission symbol rate to be equal to the sampling 
rate in the source coding part. 

The power efficiency comparison between the analog approach 
and the digital approach is shown in Fig.|4| where we deploy 10 
sensors for the digital system and we plot the power curves for 
both the case where uncoded MQAM is used (the dotted line) and 
the case where single-user capacity-achieving channel codes are 
applied to each transmitter, which gives the fundamental Shannon 
limit (the dashed line) for the digital system with orthogonal chan- 
nel usage. From the figure we see that the analog approach with 10 
nodes, which have the same observation quality and transmission 
distances as in the digital system, is more energy efficient than the 
digital approach with uncoded MQAM, but not necessarily more 
energy efficient than the fundamental limit curve. However, since 
we can support 20 nodes in the analog system with the same band- 
width requirement, with the 20 analog nodes we can achieve a bet- 
ter power efficiency (the solid line) than the digital system with 10 
nodes under the optimal single-user channel coding. This is true 
even when the extra ten nodes in the analog system have worse 
observation quality and larger transmission distances than the first 
ten nodes. 

4. CONCLUSIONS 

In this paper, we have shown that for a bounded source with un- 
known statistics and a fusion center equipped with the best linear 
unbiased estimator, we can minimize the total power consump- 
tion across all the sensor nodes under a certain distortion require- 
ment. The information collection can be implemented with an ana- 
log approach, which may be more energy efficient than the digital 
approach when only orthogonal multiple access schemes such as 
TDMA and FDMA are used. 



Analog: 20 sensors 

Digitah uncoded MQAM with 10 sensors 

- — Digital: fundamental limit with 1 O orthogonal channels 




Fig. 4. Comparison of the L^-norm of the power vector 
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